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Abstract. We present the first interferometric observations of CO(l-O) and 00(2-1) line emission from the warped 
LINER NGC 3718, obtained with the IRAM Plateau de Bure Interferometer (PdBI). This LI. 9 galaxy has a 
prominent dust lane and on kiloparsec scales, a strongly warped atomic and molecular gas disk. The molecular 
gas is closely associated with the dust lane across the nucleus and its kinematic center is consistent with the 
millimeter continuum AGN. A comparison of our interferometric mosaic data, which fully cover the ~ 9kpc 
warped disk, with a previously obtained IRAM 30 m single dish CO(l-O) map shows that the molecular gas 
distribution in the disk is heavily resolved by the PdBI map. On the nucleus the interferometric maps alone 
contain less than one half of the single dish line flux, and the overall mosaic accounts for about a tenth of the total 
molecular gas mass of ~ 2.4 x IO^A^q. After applying a short-spacing correction with the IRAM 30 m data to 
recover the missing extended flux, we flnd in total six main source components within the dust lane: one associated 
with the nucleus, four symmetrically positioned on either side at galactocentric distances of about 1.3 kpc and 
4.0 kpc from the center, and a sixth on the western side at ~ 3 kpc with only a very weak eastern counterpart. 
In the framework of a kinematic model using tilted rings, we interpret the five symmetric source components as 
locations of strong orbital crowding. We further find indications that the warp appears not only on kpc scales, but 
continues down to ~ 250 pc. Besides the sixth feature on the western side, the lower flux (a factor of ~ 2) of the 
eastern components compared to the western ones indicates an intrinsic large scale asymmetry in NGC 3718 that 
cannot be explained by the warp. Indications for a small scale asymmetry are also seen in the central 600 pc. These 
asymmetries might be evidence for a tidal interaction with a companion galaxy {large scales) and gas accretion 
onto the nucleus {small scales). Our study of NGC 3718 is part of the NUclei of GAIaxies (NUGA) project that 
aims at investigating the different processes of gas accretion onto Active Galactic Nuclei (AGN). 
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Fig. 1. Optical image (taken from the DSS survey) su- 
perimposed with the integrated CO(l-O) contours (white) 
obtained at the PdBI, with added short spacings from the 
IRAM 30 m telescope (see section IOfc lO|) . Contour lev- 
els as in Fig. ^ 

1. Introduction 

Knowledge of the distribution and kinematics of the cir- 
cumnuclear molecular gas in active galaxies is essential 
for understanding the fueling of the central engine and 
the role of gas and dust in obscuring the active nucleus 
(AGN). Although we have begun to understand the fac- 
tors that influence nuclear activity, a wealth of unan- 
swered questions still remain. In contrast to large scales 
(> 3kpc) where dynamical perturbations like galaxy col- 
lisions, mergers, or mass accretion (Heckman et al. 1986) 
as well as bars, spirals and their gravity torques (Combes 
1988) are responsible for the infall of gas, the processes re- 
sponsible for removing angular momentum at small scales 
(sub-kpc) are not very well understood. Various scenar- 
ios invoking nested bars (e.g., Shlosman et al. 1993), spi- 
rals (e.g., Martini & Pogge 1999), warped nuclear disks 
(Schinnerer et al. 2000a, b), and lopsidedness or m = 1 in- 
stabilities (Kormendy & Bender 1999; Garcfa-Burillo et al. 
2000) have been proposed. Millimeter interferometry al- 
lows qualitative studies of molecular gas emission in these 
galaxies with arcsecond spatial resolution (and even be- 
low), high spectral resolution a few kms~"'^per chan- 
nel) and high sensitivity (detectability thresholds dgas ~ 
a few tens Mqvc-^). The NUclei of GAlaxies (NUGA) 
project (Garcfa-Burillo et al. 2003a, 2003b) aims at es- 
tablishing a high-resolution and high-sensitivity CO sur- 

Send offprint requests to: M. Krips, krips@phl.uni-koeln.de 

* Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 
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Fig. 2. Positions of the IRAM 30 m observations are plot- 
ted with black crosses, and mosaic fields of the PdBI ob- 
servations are shown as circles indicating the 42" primary 
beam size for each observing field. The 30 m beamsize is 
shown in the upper left. 

vey of 12 nearby AGN covering the full sequence of ac- 
tivity types including Seyfert 1 and Seyfert 2 galaxies 
as well as Low Ionization Nuclear Emission-line Region 
(LINER) and transition objects. The CO observations are 
conducted with the IRAM Plateau de Bure Interferometer 
(PdBI: Guilloteau et al. 1992), which currently provides 
the highest sensitivity and resolution for the study of the 
distribution and dynamics of molecular gas in the nuclei of 
these galaxies. Besides case by case analyses and simula- 
tions of each object (Garcfa-Burillo et al. 2003b, Combes 
et al. 2004) , the collected data will also be used for a first- 
order approach to a statistical study of how gas flows into 
nuclei and the different mechanisms that account for fur- 
ther accretion inward. 

One of the active galaxies belonging to the NUGA sur- 
vey is NGC 3718 (Fig. 0, a peculiar, elliptical galaxy at 
a distance of 13Mpc (Pott et al. 2004; Schwarz 1985). 
Together with its companion NGC 3729 it belongs to the 
Ursa Major group. It still remains unclear if and to what 
extent a gravitational interaction is taking place between 
these two galaxies. A large warped dust lane runs across 
the entire stellar bulge of NGC 3718; it has a width of 
< 2" at the center and diverges into several smooth fil- 
aments across the bulge (1" corresponds to ^ 64 pc for 
NGC 3718). At a radius of about 1.5' from the center, the 
dust lane bends by almost 90 degrees towards the north 
and south. This warp signature is also observed in the cold 
gas distribution, with the HI line emission tracing it out 
to a radius of more than 6' (24kpc) and CO line emission 
tracing it towards the center (~ 20": Pott et al. 2004). 
The outer parts of the warp were kinematically modeled 
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2000/01 - 2001/02 - 

cont@ cont@ cont@ 

3 mm 1 mm 3 mm 

RX 11:32:34.852 11:32:34.878 11:32:34.83 

(hhmmss) ±0.003 ±0.005 ±0.01 

Dec 53:04:04.51 53:04:04.3 53:04:04.2 

(°:' :") ±0.06 ±0.1 ±0.1 

peak flux 9.4±0.7 10±2 7±1 
(mJy boam^^) 



Table 1. Intcrfcromctric positions and fluxes of the con- 
tinuum emission at 1 mm and 3 mm from the single field 
data set (2000/01) and the central field of the mosaic ob- 
servations from 2001/02 (right column). The difference in 
the 3 mm flux densities between the first epoch and the 
second is discussed in the text f section l3.1.1|) . 



calibrated relative to CRL618 and MWC349. The frequen- 
cies were centered on the redshifted CO (1-0) fine in the 
USB at 3mm and on the redshifted 1^00(2-1) line in 
the LSB at 1 mm. For each line, the total bandwidth was 
580 MHz and the spectral resolution 1.25 MHz. The inte- 
gration time for the central pointing amounts to 8hr 
on source. This gives a point source sensitivity of about 
^ TmJy 16mJy) in channels of 1.25 MHz width at 
3 mm (1mm). The precipitable water vapor ranged be- 
tween 4 and 10 mm (i.e., opacities of ~0.2-0.3) resulting 
in system temperatures of approximately 200-300 K on av- 
erage. We used the GILDAS software packages to reduce 
and map the data (Guilloteau & Lucas 2000). The synthe- 
sized beams are 4.7" x 3.5" at 41°at 3 mm and 2.6" x 2.0" 
at 45°at 1 mm. 



by Schwarz (1985), while the inner parts down to ~ 20"' 
were modelled by Pott et al. (2004). 

NGC 3718 is classified as a LINER galaxy of type 1.9 
(Ho et al. 1997). Weak, broad Ha emission with FWHM 
2350 km s~^ is emitted from the nucleus as well as strong 
[O I] A==6300Awith FWHM 570 km s^^ indicative of a hid- 
den AGN (Fihppenko et al. 1985; Ho et al. 1997). So far no 
ultraviolet emission has been detected towards NGC 3718 
(Earth et al. 1998), probably due to obscuration of the nu- 
cleus by dust, whereas a radio source has been confirmed 
at the position of the nucleus (this work; Burke & Miley 
1973; Nagar et al. 2002). 

In this paper we analyze single field and mosaic data 
obtained with the PdBI at an angular resolution of ^ 4". 
Section |21 describes the observations. In section |21 we 
present the results from the observations. A kinematic 
model is outlined in section 0] We finish with a discus- 
sion. 

2. Observations 

NGC 3718 was observed with the IRAM 30m and the 
IRAM PdBI. The IRAM 30m observations are already 
described in detail by Pott et al. (2004) and will thus 
here only be discussed in the context of the short-spacing 
correction and resolution effects. The IRAM PdBI obser- 
vations were carried out in two different modes: central 
pointing (one single field on the center) and mosaic (seven 
fields along the dust lane). 

2.1. Central pointing 

During winter 2000/2001, we mapped the CO(l-O) and 
C0(2-l) lines in a single field centered at the ra- 
dio position of the AGN in NGC 3718, i.e., Q!2ooo = 
ll'^ 32™3 4.840^ and (52ooo = 53°04'04.90" (see sec- 
tior l3.1.l1l . The PdBI was at this time deployed in the 
CD set of configurations with 5 antennae. The bandpass 
calibrator was 3C273, while the phase and amplitude cal- 
ibration were performed on 1150+497. Fluxes have been 



2.2. Mosaic 



During winter 2001/2002, a mosaic of NGC 3718 was ob- 
tained with the PdBI. Seven mosaic fields of the CO(l-O) 
emission were repeatedly observed with six antennae or a 
subset in C and D configurations. These sample the molec- 
ular gas emission along the dust lane with 12 kpc extent 
that was previously mapped with the IRAM 30 m by Pott 
et al. 2004 (Fig. 12 see also Hartwich et al. 2002 and Krips 
et al. 2003). The central field of the mosaic was set to 
the same position specified in section 12.11 and the sepa- 
ration between the pointings was 21", corresponding to 
half the primary beam size. The bandpass was calibrated 
on 3C84, 3C273, or 0923+392 with an accuracy of about 
5%. The latter source and 1150+497 were used as phase 
and amplitude calibrator, observed every 20 minutes. To 
derive the absolute fluxes of the calibrators and the target 
source, values for each observing epoch were bootstrapped 
via a standard IRAM flux monitoring program, yielding 
an accuracy of the fiux calibrations of 10% at 3 mm and 
~ 20% at 1 mm. 



The spectral resolution was set to 1.25 MHz at 
114.89 GHz, giving a contiguous bandwidth of 580 MHz. 
The integration time per field amounted to ~ 2.4 hr. The 
resulting point-source sensitivity is ^ 17mJy for a band- 
width of 1.25 MHz and 5 antennae at 3 mm. The precip- 
itable water vapor ranged between 4 and 10 mm (i.e., 
opacities of ^--^0.2-0.3) resulting in system temperatures of 
approximately 200-300 K on average. Even if the C0(2- 
1) line had been simultaneously observed in mosaic mode, 
the mosaic fields would only have been usable as single 
fields and not as complete mosaic since the separations of 
the mosaic fields are too large at 1 mm where the primary 
beam size is ^ 21". Also, the sensitivity is too poor in most 
of the fields, resulting in no detection of the C0(2-l) line 
except in the central field. 
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Fig. 3. Continuum emission of the single field observations at 3 mm (left) and 1mm (right). The cross indicates the 
position at 18 cm with EVN. Contour levels are from 2a {2a) = 1.0 (1.6) mJy bearn"^ to 9 (9.6) mJy beam~^ in steps 
of 4cr (2cr) at 3 mm (1mm). Beam sizes are indicated at lower left. 
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Fig. 4. PdBI central maps of the velocity integrated CO(l-O) (left) and C0(2-l) (right) emission in NGC 3718 (from 
velocities ^ —300 to ~ +300kms~"'^). Contour levels are from 3(7 (Sa) = 1.5 (2.7) to 5 (7.2) Jybcam^^ kms^^in steps 
of la (la) for CO(l-O) (C0(2-l)). The cross indicates the position at 18cm with the European VLBI Network (EVN; 
Krips et al. 2005, in prep.). 



3. The data 

3.1. Central emission 
3.1.1. Continuum 



Before combining the central pointing observation (single- 
field) with the new central field from the mosaic obser- 
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Fig. 5. The Spectral Energy Distribution (SED) of 
NGC 3718. IR and optical fluxes with their errors are 
shown for completeness (taken from NED). The cm fluxes 
were observed with the VLA (FIRST: Becker et al. 1995; 
Nagar et al. 2001, 2002) and thus have an angular res- 
olution comparable to the PdBI maps. The mm fluxes 
for both observing epochs are plotted (upper limit at 
1mm). The solid line indicates a fit (fjy oc i/") with 
a = -0.01 ±0.1. 



vations, the continuum for both data sets at 1mm and 
3 mm has been independently determined from the line 
free channels. In both data sets at 3 mm, i.e. single field 
and central mosaic field, we detect faint continuum emis- 
sion in the central region of NGC 3718. At 1mm, only 
the single-field data set reveals continuum emission (see 
Fig. O; we do not find any emission in the central mo- 
saic field probably due to sensitivity limitations for the 
1 mm mosaic data. At 3 mm as well as at 1 mm, the contin- 
uum emission is point-like relative to the respective beams 
(compare Table The continuum declines by ^ 30% be- 
tween the two observing epochs at 3 mm (see Table 
Assuming an upper limit of 3a = 7.5 mJy beam~^ at 1 mm 
for the continuum of the central field from the mosaic, we 
obtain a decline of > 20% relative to the single field data 
(flux lOmJybeam"^) consistent with the 3 mm data. 
As the atmospheric conditions were quite similar (similar 
amounts of water vapor) between both data sets, we can 
mostly exclude the possibility of atmospheric absorption 
which could cause an artificial drop of the fiux between 
the data. The phase and amplitude calibrators were never 
separated from NGC3718 by more than a few degrees in 
elevation. The ^ 10% uncertainty in 3 mm fiux calibration 
cannot explain all of the observed variability, although this 
is not the case at 1 mm. We do not find any support for 
a variable continuum source at other wavelengths in re- 
cent EVN and MERLIN measurements at 18 cm and 6 cm 
(Krips et al. in prep., epochs 2001 to 2004) or other cm 
observations (at 2cm and 3.6cm; e.g., Nagar et al. 2002). 
In conclusion, the detected mm-variability might be an 
intrinsic property of NGC 3718 but has to be further in- 
vestigated. 

The derived positions are consistent between the two 
observing epochs as well as between 1 mm and 3 mm. They 
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Fig. 6. CO(l-O) {upper panel) and C0(2-l) {upper panel) 
spectra of the central emission in NGC 3718 observed 
with the PdBI (dark grey; integrated over the emission 
area), the IRAM 30m telescope (light grey; Pott et al. 
2004) . The short-spacing corrected spectrum is also shown 
(dashed line; integrated over the emission area). Spectra 
are continuum subtracted. The rms noise level is here 
- 4mJy (- 6mJy) for the PdBI data for CO(l-O) (C0(2- 
1)) (Aw=13kms"ifor both). 



also agree with the positions of the nucleus obtained at cm 
and optical wavelengths within the errors. Accounting for 
a weak variability, a comparison between our mm fluxes 
and cm fluxes observed with the VLA in 1997 (Becker et 
al. 1995), 1998 (Nagar et al. 2001), and 1999 (Nagar et 
al. 2002) indicates a spectral index of -O.OliO.l depend- 
ing on which of our two observing epochs is used. This 
corresponds to a flat to inverted spectrum for the central 
source in NGC 3718, like those found for many other low- 
luminosity active nuclei in Seyfert galaxies and LINERs 
(Wrobel & Heeschen 1991; Slee et al. 1994; Nagar et al. 
2000). The flat radio spectra exclude emission from op- 
tically thin synchrotron emission, as is often seen from 
extended radio jets, since this produces steep radio spec- 
tra. Thus, synchrotron self-absorption or advection dom- 
inated accretion flows (ADAFs: Narayan, Mahadevan & 
Quataert 1998; Quataert 2001) have to be considered. A 
compact jet in combination with an ADAF (Falcke 1996; 
Falcke & Marloff 2000; Falcke 2001) also cannot be ex- 
cluded. If a longitudinal pressure gradient were to cause 
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the jet to accelerate rapidly, the integrated radio emission 
of the compact core jet would also have a slightly inverted 
radio spectrum. 

3.1.2. Line 

After subtracting the continuum emission in the uv plane 
from the (continuum+line) channels, CO(l-O) and C0(2- 
1) line emission is detected in the central ^ 10" (= 1 kpc) 
of NGC 3718. The combined central maps are shown in 
Fig- El In both lines the emission is slightly extended. The 
positions of the integrated CO(l-O) and C0(2-l) emis- 
sion are similar to the nuclear position derived from the 
radio. Also, the iso- velocity diagram in Fig. shows that 
the dynamical center of the CO emission lies on the radio 
nucleus. A comparison of the flux observed with the PdBI 
and the IRAM 30 m single dish reveals that the interfer- 
ometer recovers less than one half of the peak line flux 
seen with the 30 m telescope (see Fig. . This indicates 
that the emission in NGC 3718 must be much more ex- 
tended and diffuse than is apparent in the integrated map 
of Fig. m Assuming a standard Mgas(=M(H2-hHe)=1.36 
xM(H2)) to ico(i-o) '^atio of 5A^0(Kkms-ipc2)-i 
(Solomon & Barrett 1991), the central molecular gas mass 
Mgas amounts to 1 x IO^Mqih the PdBI map, in con- 
trast to 3 X lO^TWoin the 30 m map (Pott et al. 2004). 
Downes et al. (1993) emphasize, however, that the stan- 
dard galactic conversion factor between gas mass and CO 
luminosity might fail for the centers of galaxies where the 
molecular gas is bound by the total gravitational poten- 
tial of the galaxy rather than by self-gravity. In such a 
case, the conversion factor would be up to 5 times lower 
than the standard galactic one, resulting in gas masses of 
~ 0.2 X lO'^Xofor the PdBI map and 0.6 x lO'^TWofor 
the 30 m map. 

By comparing the distributions of the central C0(1- 
0) and C0(2-l) line emission, we find that the C0(1- 
0) emission shows an elongation towards the south, al- 
most perpendicular to the beam. We estimate the inte- 
grated line ratio R21 =CO(2-1)/CO(1-0) to - 0.6 indi- 
cating subthermal excitation conditions and cold gas. To 
obtain this value, we have first truncated the inner uv re- 
gion of the CO(l-O) data and the outer uv region of the 
CO (2-1) data to get comparable uv coverages and resolu- 
tions for the two transitions, integrated both lines over the 
same velocity range of -300kms~^to -|-300kms~^and then 
transformed the fiux (Jansky) into temperature (Kelvin). 
However, due to the limited signal-to-noise and the differ- 
ence of the weighting factors between both transitions, the 
reliability of this estimate is still questionable and thus the 
value of ~ 0.75 derived in the center by Pott et al. (2004) 
might be a better indicator for the physical conditions in 
NGC 3718. 

Fig.Qshows the zero order moment maps of both tran- 
sitions overlaid with contours from the first order moment 
maps and the position velocity cuts taken along the kine- 
matic major axis. The CO (2-1) maps are quite noisy but 
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Fig. 10. PdBI mosaic map of the velocity integrated 
CO(l-O) emission in NGC3718 (upper figure; from -300 
to -|-300kms~^). Spectra are shown in the lower panel. 
Contour levels are from 2a =1.0 Jybeam^^ knis^^to 
4 Jybeam^^ kms^^in steps of la. 
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Fig. 7. Mean velocity fields of the CO(l-O) {upper left) and C0(2-l) {lower left) lines with their respective position- 
velocity maps (CO(l-O): upper right; C0(2-l) lower right) taken along the axis indicated in the maps {solid line: major 
axis). The grey scale in the images on the right is the integrated line emission of CO(l-O) and C0(2-l) respectively. 
The light grey stars in the lower right panel illustrate the positions of the CO(l-O) components identified in the upper 
right panel. Left: Velocity contours are in steps of 25kms~^for CO(l-O) and for C0(2-l). Right: Contour levels are in 
steps of 10% from the peak starting at 30%. The position velocity maps were Hanning smoothed. 



are plotted for consistency and comparison purposes. The 
position velocity diagram of the CO (1-0) line emission re- 
veals a quite extended {~ 8") velocity gradient between 
±150kms^^. Two further features appear at very high 
velocities (~ ±220kms~^) with a somewhat steeper gra- 
dient (~ ±2"). Consistent with Pott et al. (2004), the 
global velocity range of the inner ^ 15" amounts to al- 



most ±220kms~'^. However, this "central" feature in the 
position velocity diagram extends in Pott et al. (2004) out 
to a radius of 15-20", while in our interferometric data it 
stops at a radius of < 10". This rapidly rotating nuclear 
component will be discussed in section 1^31 The line emis- 
sion shows a strong asymmetry in the position velocity 
diagram with only little emission towards the blueshifted 
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Fig. 8. PdBI central maps corrected with 30m data for short-spacing effects of the integrated CO(1~0) (left) 
and C0(2-l) (right) emission in NGC 3718 at same size. Contour levels are from Str (Str) = 2.0 (2.4) to 8.0 
(10.4) Jy beam~^ kms~"'^in steps of Str {2a) for CO(l-O) (C0(2-l)). The white cross indicates the position of the 
AGN radio continuum source at 18 cm derived from EVN data. 
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Fig. 9. Comparison of pv-diagram along axis shown in Fig. [7| between PdBI-only data for the central pointing (left) 
and PdBI-|-30m data for the central mosaic {right). The asymmetry is lowered in the central mosaic map but does not 
disappear totally. Thus, it cannot be totally due to missing short spacings. Contour levels are in steps of 10% from 
the peak starting at 30%. The position velocity maps were Hanning smoothed. 



velocities. Section will show that this asymmetry can- 
not be totally traced back to missing short spacings. 
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Component 


R.A. (offset) 


Dec. (offset) 


Flux 




["] 


["] 


Jykms"^ 


SE2 


+52±3 


-34±3 


1.2 


SEl 


+20±3 


-8±3 


1.3 


C 


Oil 


Oil 


3.4 


NWl 


-19±1 


+9±1 


2.9 


NW2 


-49±1 


+36±1 


2.7 



Table 2. Positions and fluxes for CO(l-O) components 
in the PdBI+30m map. Offsets are with respect to 
a20oo=llh32m34.8 and (52ooo=53°04'04.9". 



3.2. Overall emission 

3.3. PdBI-only data 

Outside the central mosaic field, there is no evidence of 
continuum emission in our data. The 30 m map of Pott et 
al. (2004) has already revealed the existence of two addi- 
tional maxima in the CO (1-0) emission which are offset 
by ~ 30 — 40" from the nucleus. We find further maxima 
in the new PdBI mosaic map which are located farther 
out at ~ ±60" (compare Figs. Cni and HH). In Fig. ^ 
we have labeled the eastern components SEl and SE2, 
the western components NWl and NW2, and the cen- 
tral feature C. There are also indications for further max- 
ima between NWl and NW2 and also between SEl and 
SE2 but they still lack enough sensitivity (< 2a) in the 
PdBI-only map to be further substantiated at this point. 
The brightest line in the spectra turns out to be NW2 
and C, whereas the other three features appear to be 
rather weak with a signal-to-noise ratio at line peak of 
^ 4 (Fig. I10|l . However, NW2 appears to be quite weak 
in the integrated map which might be traced back to its 
weak line width of only '^50kms~^. The integrated map 
was produced by averaging the channel maps over a total 
velocity width of ~600kms~^. A comparison between the 
global fluxes derived from the 30 m and the PdBI data 
again shows that a large part of the - diffuse, extended - 
flux (at least ~ 70% of the total flux) must be resolved by 
the interferometer due to missing short spacings. Adding 
up all the PdBI spectra leads to a total molecular gas 
mass of A/gas ^ 2.5 x 10'^ Mq. This is a tenth^ of the to- 
tal molecular gas mass derived from the 30 m spectra of 
~ 2.4 X 10^ A40(Pott et al. 2004). 

3.4. PdBI+30m data: short spacing correction 

To restore missing extended and diffuse emission, we com- 
bined the CO(l-O) data from the IRAM 30m telescope 
with our interferometric maps. We varied the weights at- 
tached to the 30m and PdBI data to flnd the best com- 
promise between good angular resolution and almost com- 
plete restoration of the missing extended flux. In the 
end, a factor of 4 was applied to the weights of the 30m 
data, leading to the recovery of ^ 90% of the missing 



^ Besides the decreased flux measured with the PdBI also 
the line width is lower by < 20% with respect to the 30m data. 
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Fig. 12. Tilt angles uj{t). The old model of Pott et al. 
(2004) is given as a solid line, while a dashed line is used 
for our new model. 
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Fig. 13. The integrated intensity map of our model is plot- 
ted here, before spatial filtering to match the response of 
the interferometer and single dish telescope. 

flux (see below) at a still reasonable angular resolution 
of 7.8x4. 5"@20°. A higher weighting factor would have 
increased the image restoration to over 90% but at the 
same time also decreased the resolution to ~ 10", while 
a lower weighting factor would have had a higher an- 
gular resolution but a worse image restoration due to a 
lower signal-to-noise ratio. The combined maps of 30m 
and PdBI data have finally been cleaned using standard 
program packages provided by the GILDAS software for 
mosaic data. The results are presented in Fig. ISl fell II for 
the single field and mosaic maps. In Fig.|2| we compare the 
position- velocity diagrams along the axis used in Fig. |7|for 
the PdBI-only data and the PdBI-|-30m data. The asym- 
metry discussed in Section 13.1.21 is still present even if 



slightly weakened and hence cannot be totally traced back 
to missing short spacings. This indicates that the asymme- 
try seems to be an intrinsic property of NGC 3718. A first 
comparison between the PdBI-only and the PdBI-|-30m 
maps shows in the mosaic maps that the short spacing 
correction strengthens the inner components (NWl, C 
and SEl) that appear to contain most of the extended 
and diffuse emission. If we now derive the total molecu- 
lar gas mass, we obtain ~2.1xlO*A^o, which is ^^90% 
of the value found for the 30m-only data. Besides the 
five main emission peaks seen in Fig. ^1 further emis- 
sion between NWl and NW2 and SEl and SE2, as al- 
ready indicated in the PdBI-only map (Fig. is sub- 
stantiated in the PdBI-|-30m map. Their gas masses are 
estimated to ~ 0.2 x 10* A^©, respectively. Globally, the 
western part of the CO (1-0) emission appears to contain 
more gas (Mgas — 10* A^Qwithout C) than the eastern one 
(Mgas ^ 0.5 X 10* A^Qwithout C). The CO(l-O) emission 
cannot be associated with a single straight line but rather 
delineates an "S-like" shape. 

The next sections will deal with the modeling of the 
complex gas distribution in NGC 3718 for a better under- 
standing and a more quantitative interpretation. We will 
also examine in detail the effects of missing short spac- 
ings on our data by applying our procedure to a model 
and comparing the results thus obtained with the obser- 
vations. 

4. Modeling the CO distribution 

The special characteristics of the HI data for NGC 3718 
led to the conclusion (Schwarz 1985) that the atomic gas is 
warped at radii > 8 kpc in a configuration similar to that 
of the dust lane. IRAM 30 m observations of the molecular 
gas by Pott et al. (2004) showed this warp extending down 
to smaller radii of ~ 1 kpc. Previous modelling based on 
the Pott et al. CO and the H I data can hence be fine-tuned 
with the new information gained from our significantly 
increased spatial resolution- ^4 — 7" for our PdBI-only 
map vs. 21" for the Pott et al. (2004) 30 m map. 

4.1. The model 

Our simulations of the mosaic maps are based on a tilted 
ring model similar to that of Pott et al. (2004); that pa- 
per presents more modelling details. Since the Pott et 
al. (2004) model of NGC 3718 reproduces the 30 m ob- 
servations with only 3 main maxima, we had to modify 
the parameters to adapt their model to our PdBI results. 
The revised best-fit parameters were obtained by start- 
ing with the model of Pott et al. (2004), propagating it 
through the PdBI+30m response function (described in 
section l^^ . comparing the result with the observed data 
(PdBI-|-30m), and finally tuning the parameters until the 
match is optimal. In the new best-fit model, we mainly 
use different tilt angles w(r) (of a given ring relative to 
the assumed plane of the innermost reference disk) for the 
inner radii (< 50"). Our inner rings are less warped, i.e., 
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Fig. 14. Model vs. observation. Integrated maps of the 
PdBI data only (upper panel), the PdBI and 30m data 
combined (middle panel) and the 30 m data only [lower 
panel). The model (grey scale) is overlaid with solid con- 
tours of the observed data. Contour levels: PdBI-only 
2(7=1.0 to 4Jybeam-ikms-iby la; PdBI+30m 2cr=1.0 
to 6.0 Jy beam~^ kms~^by Icr; 30m-only 3(t =1.0 to 10 
by 1.0 Jybeam"^ kms~^. 



the tilt angles are smaller for radii from 40" to 60" than 
in the previous model, but they remain constant down to 
radii of 1" (i.e., in our model, the warp continues down 
to a radius of ~ 1"; see Fig. ll2f) . Furthermore, we use a dif- 
ferent intensity profile from Pott et al. (2004). The earlier 
work assumed an integrated line flux distribution for the 
rings that is very bright at the center but then decreases 
rapidly, while the intensities in our model show a slower 
decrease with radius. We also adopt a different inclina- 
tion (70°instead of 60°) and position angle (-60°instead 
of -85°) for the galactic disc in space. The rotation curve, 
however, remains the same for the new model. The inte- 
grated map of our model is shown in Fig. E| (the model 
has not yet been filtered through the interferometer or 
single dish response function). While the central (C) and 
the two outer (SE2, NW2) components are quite compact 
and bright, the two inner features (SEl, NWl) are rather 
weak and relatively extended on ^ 10" scales. 



4.2. Simulating the short-spacing correction 

To allow for a direct comparison between the derived 
model and the CO data (PdBI-only, PdBI-hSOm, 30m- 
only), we simulated the effect of the different observations 
on the model. To determine how our model galaxy would 
look when observed with the interferometer (before and 
after short-spacing correction), we first produced the sim- 
ulated PdBI map with a program that computes an inter- 
ferometric image (upper panel of Fig. ll4|l from a given uv 
coverage (taken here from the observations) and the given 
three-dimensional intensity distribution (taken here from 
our best-fit model shown in Fig. I13|l. Only the most com- 
pact features are visible in the PdBI-only map while the 
extended emission is resolved, resulting in a loss of around 
60% of the original flux. These compact features corre- 
spond to the compact orbital crowding points SEl, SE2, 
NWl and NW2 ( Fig. I13|l . The remaining emission from 
the tilted rings appears to be extended and mainly con- 
centrated around SEl and NWl. It cannot be seen with 
the PdBI alone due to spatial filtering, meaning that only 
the five compact crowding points remain in the PdBI-only 
map. The extended emission around SEl and NWl is thus 
heavily resolved (compare Fig. ll3|l . The best-fit model was 
also smoothed to the beam (~ 21") of the 30 m telescope 
and the spectra were taken at the same positions as for 
the observations. Accounting for normalisation, these data 
were transformed into the 30 m file that would have been 
produced if such a distribution had been observed with 
the IRAM 30 m telescope. This file is thus suitable for 
comparison with the measured 30 m data (see lower panel 
in Fig. I14|l . We indeed find only three main maxima, a 
merging of the inner two (SEl, NWl) with the outer two 
maxima (SE2,NW2) and an S-line morphology. 

As a final test, we combined the simulated 30 m and 
PdBI maps. All three simulated maps are plotted in 
Fig.EKffre?/ scale). The loss of flux due to the reduced sen- 
sitivity to diffuse extended emission in the interferometer- 
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are in steps of 10% from the peak starting at 20%. The position velocity maps were Hanning smoothed. 



only data is indeed reproduced by the model. A com- 
parison of the simulated peak fluxes from the PdBI, the 
PdBI+30m, and the 30 m data with each other suggests 
that the PdBI has lost a factor ~ 2 — 3 of the peak line 
flux, consistent with what is observed. The short-spacing 
correction not only reproduces the correct fluxes but also 
the correct distribution without any artefact or inconsis- 
tency in excess of 1% for the model. 



4.3. Model vs. Data 

We have carried out the same steps to combine the ob- 
served PdBI and 30 m data fsection l3.4|l . The results are 
plotted in Fig. 1141 as contours over the simulated maps. 
We find a very similar behavior in the observed and model 
maps. Besides the inner two maxima which are more easily 
visible, the flux level also increases. Thus, a better distinc- 
tion between signal and noise is possible. As the two east- 
ern components are less dominant (by a factor 2) than the 
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western components in the integrated map over the whole 
velocity range, we show separate plots for each compo- 
nent integrated over the appropriate velocity range (see 
Fig. I15|l . Such an asymmetry of the gas distribution has 
already been found in many other galaxies as well (e.g., 
Richter & Sancisi 1994) and cannot be explained by our 
symmetric tilted ring model. It might be an indication for 
a tidal interaction with the companion galaxy which would 
result in an asymmetric gas distribution. Thus, our warp 
model could presumably be improved by adding asymme- 
tries to account for the putative interaction. 

The maps of the data [left side) and the model [right 
side) are separately shown in Fig. 1161 along with the good 
agreement between the model and the data in terms of 
spectra taken at each of the five components. A somewhat 
weaker eastern component is also visible in the model. 
This is due to the asymmetric sampling of the 30 m- 
spectra with respect to the CO centroid. Simulated and 
observed line widths are consistent within the noise. In 
Fig. El we have derived first order moment maps with 
a 3(7 clipping over a velocity range of ~ — 300km s~^to 
-|-300kms^^and taken position- velocity cuts along the ma- 
jor axis of the galaxy for the observed and simulated data. 
Compared to the pv-diagram of a standard gas disk, a 
more complex structure is seen in the pv-diagrams of the 
observations as well as of the simulation due to orbital 
crowding effects. As might be expected for a rapidly rotat- 
ing nuclear disk with a size of 1.5 kpc, a steep linear rise 
is indicated in the position-velocity diagram. This feature 
agrees well in orientation and size with the central field 
velocity component visible in the CO(l-O) map of Fig. [3 
We again find that this central feature in the position- 
velocity map extends only over a range of ±10" in Fig. 1171 
which is smaller than the interval derived by Pott et al. 
(2004). The difference from the position velocity diagram 
in Pott et al. (2004) has a logical and simple explanation. 
In their maps, {SEl, SE2} and {NWl, NW2} merge (in 
the 30 m beam) into almost single SE and NW features, 
respectively, and thus all lie on the major axis. In our in- 
terferometric map, SEl and NWl are no longer aligned 
with this major axis and thus appear only weakly in the 
position- velocity diagram of Fig. El The high velocity gas 
returns if the major axis is rotated by about '--^10°. Wc thus 
have three different velocity regimes. The first ranges from 
±10" with Av ~ ± ~ ISOkms^^and can be identified 
with the gas occupying the inner rings. The second region 
extends to ±30" with Aw ±220kms~^and is caused 
by intermediate rings. The third extends to ±70" with 
Av = ±220kms~^and originates from the outer rings. 

5. Discussion & Conclusions 

The paper presents an analysis based on arcsec angular 
resolution PdBI observations of the gas distribution in 
NGC 3718, complementing previous studies conducted by 
Pott et al. (2004: CO on > 20" scales) and Schwarz (1985: 
HI on > 30" scales). After making a short-spacing correc- 
tion using the IRAM 30 m data (Pott et al. 2004), we have 



modelled and interpreted the mosaic observations of the 
total emission as well as single field observations of the 
central kpc. Three compact and three extended compo- 
nents were found in the mosaic map of the CO(1~0) emis- 
sion. At the higher angular resolution of the PdBI, the gas 
structure turns out to be more complex than previously 
appreciated. To explain the detected molecular gas distri- 
bution, we were forced to modify the Pott et al. tilted ring 
model slightly. The revised model now not only reproduces 
well the PdBI+30m observations, but due to the higher 
resolution of the observations, also places much tighter 
constraints on the warp. The disk is most likely already 
warped on arcsec scales (~ 4" = 250 pc) and is evident 
within a radius of 20" (the S-like shape between SEl, C, 
and NWl) even if the warp is less dominant within 50" rel- 
ative to the original Pott et al. model. The five symmetric 
CO(l-O) peaks (SE2, SEl, C, NWl, and NW2) detected 
by the PdBI can be explained by orbit crowding effects. 
The two outer maxima merge into one feature on each 
side of the center in the beam of the IRAM 30 m telescope 
and thus at this resolution appear as single components. 
Besides the gas distribution, the gas kinematics can be 
almost completely explained by the model. The position- 
velocity diagrams of both data sets - mosaic and central 
pointing - unveil rapidly rotating nuclear gas within a ra- 
dius of ^ 700 pc, as seen in the simulated data. The high 
velocities in the center are produced by the inner warped 
rings. Thus, the dynamics of the molecular gas can be 
traced back to a continuous warped disk with concentric 
orbits. 

As in Pott et al. (2004), our model does not require 
a bar to explain the observed gas properties down to 
~ 250 pc. However a warp alone is not sufficient to ac- 
count for all data, it must be combined with asymmetry 
in gas density between the eastern and western side. This 
asymmetry is logically coming from the same mechanism 
producing the warp: i.e. an interaction with the companion 
galaxy. The asymmetric gas distribution within the cen- 
tral 250 pc provides the possibility for an accretion of gas 
onto the central engine. However, the transport of gas to 
the black hole still raises many questions and problems. 
Even in the Milky Way, where a warped circumnuclear 
molecular torus (Guesten et al. 1987) provides a reservoir 
for accretion at a distance of only ^ 1 pc from the nu- 
cleus, the situation is not yet understood. In this context, 
the mechanisms of gas fueling to the inner few pc/subpc 
in active nuclei appear even more complex. Besides the 
comprehension of our own Galactic center, only detailed 
studies of further nuclei with the highest angular resolu- 
tion can provide answers to the wealth of open questions. 
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